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Abstract 



This article describes improved measurements by CLEO of the DfD* 
and D*~^D*~ branching fractions, and first evidence for the de- 

cay 5+ Di*^^D**^, where L>**° represents the sum of the Z)i(2420)°, 
L>|(2460)'', and Di{j = 1/2)'' L = 1 charm meson states. Also reported is the 
first measurement of the D*'^ polarization in the decay 5° A par- 

tial reconstruction technique, employing only the fully reconstructed and 
slow pion TT^ from the D*~ D^t^J decay, enhances sensitivity. The observed 
branching fractions arc B{B^ D+D*-) = (1.10 ± 0.18 ± 0.10 ± 0.28)% , 
^ D*+D*-) = (1.82±0.37±0.24±0.46)% , and B{B+ oi*^'^ D**°) = 
(2.73 ±0.78 ±0.48 ±0.68)%, where the first error is statistical, the second sys- 
tematic, and the third is the uncertainty in the Df — >■ ^tt"*" branching fraction. 
The measured longitudinal polarization, F^/r = (50.6 ± 13.9 ±3.6)%, is 
consistent with the factorization prediction of 54%. 
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Typeset using REVT^ 



I. INTRODUCTION 



Measurements of weak decays of B mesons are fundamental to testing and understanding 
the standard model. Previous measurements of the inclusive B — »• DfX branching fraction 
report a value of (12.1±1.0±3.0)%. The first error is the combined statistical and systematic 
uncertainties, and the second is due to the uncertainty in the — > ^vr"*" branching fraction. 
This is significantly larger than the sum of production from exclusive b — > ccs modes 
observed to date These exclusive modes, of the form B D^D, B D*~^D, B 
DfD*, and B D*+D\ sum to (6.6 ± 1.3 ± 1.7)% for the S+ case and (4.8 ± 1.0 ± 1.2)% 
for the 5°. This yields a deficit of (5.5 ± 1.6)% for the S+ and (7.3 ± 1.4)% for the 
5°, where the 07r^ branching fraction uncertainty does not affect this difference 

This article reports new measurements of i? — > D(*)+_D*(*) decays from CLEO.Q First 
evidence is offered for the decay B^ ]ji*)+f)**o^ where D**^ denotes the sum of the 
Z)i(2420)°, 1)^(2460)°, and Di{j = 1/2)° L = l charm meson states. This decay mode may 
bridge a substantial portion of the inclusive and exclusive rate difference. Also reported are 
improved measurements of the modes B^ —>■ DfD*~ and B^ D*^D*^. These decays 
occur predominantly via the spectator diagram of Figure |1]; the decays into a or 
D*~^ meson, and the charm anti-quark and spectator quark hadronize as either a D* or D** 
meson. 
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FIG. 1. The spectator diagram for B decay. 



Additionally, this Article presents the first measurement of D*^ polarization for the mode 
B^ — s> D*^D*^, providing an effective test of the factorization assumption in S — D*^X 
decays with high g^, where = Mx, and X is a vector meson. Factorization assumes 
the lack of final state interactions between the products of hadronic B decays, and has 
successfully predicted the vector- vector polarization of the low mode B D*~ p It 
is possible that the factorization assumption of no final state interactions may be simplistic 



^Reference to a specific state or decay includes the charge-conjugate state or decayThe notation 
D^s^~^ in this context means either a or a denotes the sum of D* and D**, and D** 

denotes the sum of the charged and neutral D**^ states, the specifics of which are discussed 

in Section |IV A|. In shortened form, DgD* denotes DfD*^, D*D* denotes D*~^D* 

■)(*)+^**0 



denotes the sum of Di*^^D**- and D 



and Di*^D** 
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and inapplicable to modes of higher such as D*~^D* ; however, the results presented 

here are consistent with the factorization prediction. 

Previous measurements of 5° ^ DfD*- and 5° D*+D*- at CLEO and ARGUS 
made use of the full reconstruction technique |l|] 0, which requires reconstruction of all 
particles in the final state. The most recent CLEO results using full reconstruction reported 
relatively small event yields of 18.4±4.5 and 17.7±4.4 in the D^D*~ and D*^D*~ channels, 
respectively. Following these, a partial reconstruction technique was developed that required 
only some of the —>■ D*~^D*~ final state particles, reporting an increased sample size of 
76± 11 events §. 

This analysis employs a more refined partial reconstruction technique, using only the D'^ 
and the soft pion vr^ from the D*^ D^t^^ decay, thereby increasing the statistics over 
full reconstruction by a factor between five and eight, depending on mode. The analysis is 
sensitive to any B D+D*-X final state, such as B ^ D'f^+D**, when D** D*-7i. The 
method is based on techniques developed by CLEO for improved measurement of — > ^vr"'" 
@ and 5 ^ D*TT g. 

After a short description of the detector and the criteria used for selecting charged 
particle candidates in Section the and tt~ reconstruction is described in Section 



In Section ^ the partial reconstruction technique is developed for separating the combined 
£){*)£)*{*) signal from background. Once the background levels have been determined, in 
Section a two-dimensional parameter space is defined and used to separate the individual 
DsD*, D*D*, and signals, followed by a review of systematic errors in Section [V^. 



The polarization of D*D* production is measured and compared with the factorization 
prediction in Section |V11| , and the results summarized and discussed in the final section. 



II. EVENT SELECTION 

The data used in this analysis were collected at the Cornell Electron Storage Ring (CESR) 
between 1990 and 1995, and consist of hadronic events produced in e~^e~ annihilations. The 
integrated luminosity of this data sample is 3.14±0.06 fb~^ collected at the T(4S') resonance 
(referred to as on-resonance data), and 1.69 ± 0.03 fb~^ from a center-of-mass energy just 
below the threshold for producing BB mesons (referred to as off-resonance or continuum 
data). The on-resonance data corresponds to (3.36 ± 0.06) x 10^ BB pairs. 

The CLEO 11 detector is used to measure both neutral and charged particles with ex- 
cellent resolution and efficiency []TU|. Hadronic events are selected by requiring a minimum 
of three charged tracks, a total visible energy greater than 15% of the center-of-mass en- 
ergy (this reduces contamination from two-photon interactions and beam-gas events), and 
a primary vertex within ±5 cm in the z direction and ±2 cm in the r-0 plane of the beam 
centroid. 

Charged tracks are required to be of good quality and consistent with the primary vertex 
in both the r-0 and r-z planes. Tracks must also have dE/dx and time-of-fiight information 
consistent with their pion or kaon hypotheses, when such information exists and is of good 
quality. 

Apart from the visible energy criterion, neutral particles were not used in this analysis. 



A GEANT |TT| based Monte Carlo simulation was used to generate large samples of the 
individual signal modes from T(4S') BB decays, and model their interactions 
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with the CLEO detector. These samples were then processed in the same manner as the 
data. Further discussion of the simulation is given in the treatment of systematic errors. 

III. AND SLOW TTg RECONSTRUCTION 

The is reconstructed through the 07r+, K^K^ decay channel, which has 

a signal-to-background ratio nearly two times higher than the next cleanest decay mode 
0. Fast Tx^ /K^ tracks {p > 200 MeV/c) must originate within ±5 cm in the z direction 
and ±5 mm in the r-0 plane of the beam centroid. For slow Tr"*"/ tracks {p < 200 MeV/c) 
the z requirement is loosened to within ±20 cm. The K^K~ invariant mass is required to 
be within 9 MeV of the mass. Two angles are used in suppressing background. The first is 
the decay angle 6^, which is the angle between the (p direction in the D'^ rest frame and 
the boost direction. Requiring cosOd < 0.80 eliminates a large combinatoric background 
peak near cos On = 1 resulting from the numerous low momentum pions, while the signal 
is constant in cos 9d- The second angle is 9h, the decay angle between the and D'^ 
direction in the rest frame. Due to the helicity the signal follows a cos^ 6h distribution 
while the background is constant in cos 6h- Requiring | cos6h\ > 0.35 removes 35% of the 
background and retains 96% of signal. The resulting 07?"*" invariant mass spectrum is shown 
in Figure ^ and the 07r+ mass is then required to be within 12 MeV of the D'^ mass. 
Finally, the kinematics of -B ^ decays constrain the magnitude of momentum 

to between 1250 MeV/c and 1925 MeV/c, and these requirements are imposed here. 
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FIG. 2. The (j)TT~^ mass spectrum for the on-resonance data. The 07r^ mass is further required 
to be within 12 MeV of the D+ mass. 
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The slow pion vr" from the D*~ must have charge opposite to the and originate 
within ±5 mm in the r-0 plane of the primary vertex. No z requirement is placed on the 
7r7, but it must have a momentum greater than 50 MeV/c and less than 210 MeV/c. 

IV. SEPARATION OF SIGNAL FROM BACKGROUND 

A. Two-Body B decays to DfD*- Final States 

At the CLEO II experiment, e"'"e~ collisions can create an T{4S) resonance, which decays 
to a pair of B mesons. The -B's are produced nearly at rest (/3 = 0.0646) and, for the decay 
chain T(4^) ^ 5° ^ D+D*-, and D*- DH;, the Df and soft pion vr," are 

nearly back-to-back in the lab frame because of the small 5.83 ± 0.03 MeV energy release in 
the D*^ D^T^^ transition. By making use of their relative direction, as well as the beam 
energy and kinematic constraints of the decay, the and the vr" allow reconstruction of 
the DfD*~ final state. 

Other two-body B decays leading to DfD*^ final states, with strong (Z}+,7r^) correla- 
tions, are summarized in Table |. These are modes producing a D*~^ that decays to Df'j 
or Df-K^, or producing a D** that decays to D*^7t. It should be noted that this method is 
not sensitive to -B ^ D**D*^, as the D** decays predominantly to DK and no D** DgX 



decays have been observed [O]. Other relevant modes, such as three-body B decays of the 



form B Di*^+D*- 'TT, are treated in the discussion of systematic errors. 

TABLE I. DfD*~^ final states from two-body B decays. 



B Decays B^ Decays 

B° D*+D*- 
(where D*^ — > 
Z)+7 / D+n^) 

B^ D+D**- B+ D+D**^ 

B^ D*+D**- B+ D*+D**° 

and D*+ and D*+ 



B. D** Properties and B Dy'^D** Decays 

In this measurement several different D** states contribute to the B [)*(*) decays. 

The relevant D** characteristics are summarized here, beginning with the neutral D**^ 
which, as an L = 1 charm meson, represents four distinct quantum states. Two of these 
states, the Di(2420)° and 1^2(2420)°, have been characterized by experiment as relatively 
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narrow resonances 



mi]. The two other states, the Di{j = 1/2)° and Do{j = 1/2)°, are 



expected to be much broader 0]. A prehminary first observation of the Di{j = 1/2)°, 
confirming its broadness, was recently reported by CLEO [|15|. Although the -Do(j = 1/2)° 
remains experimentally undetected, conservation of parity and angular momentum forbids 
its decay to D*^7i^, so it does not contribute to this measurement. Table |I| gives the masses, 
widths, and allowed decays of the neutral D*' 



TABLE II. D**^ Properties 



State 




Mass (MeV) 


Width (MeV) 


Allowed 
Decays 


D*oU = 1/2)° 


0+ 


Not Yet Observed 




Dtt 


L>i(2420)° 


1+ 


2422.0 ± 1.8 


18-9+31 


D*TT 


Di{j = l/2r 


1+ 


246lli 


2901^°^ 


D*TT 


£'1(2460)0 


2+ 


2458.9 ± 2.0 


23 ±5 


DTT,D*Tr 



In accordance with current experimental limits, the masses and decay widths of the 
charged Di{j = 1/2) ~, i5i(2420)~, and D2(2460)~ are assumed identical to their corre- 
sponding neutral D**^ counterparts [|l^. Like the -Do(j = 1/2)°, the -Do(j = l/2)~ does not 
decay to D*~'Tr. 

Throughout this Article D**+ denotes the sum of the charged L'i(2420)+, £'^(2460)+, 
and Di{j = 1/2)+ states, while D**° denotes the sum of the neutral Di(2420)°, D^(2460)°, 
and Di{j = 1/2)°, and D** denotes the sum of the three D**~^ and three D**° states. 

Conservation of isospin and angular momentum predicts the branching fractions for the 
J = 1 charged and neutral D** D*~n decays. Heavy quark effective chiral perturbation 
theory evaluates the branching fractions for the J = 2 case [§. 



B{D^{j = 1/2)- - 


D* 


-7r°) = 


= 1/3, 


(1) 


^(^1(2420)- - 


D* 


-7r°) = 


= 1/3, 


(2) 


B{D;{2m)- - 


D* 


-7r°) = 


= 1/10, 


(3) 


B{D,{j = l/2f- 


D*- 


-7r+) = 


= 2/3, 


(4) 


i3(Di(2420)° - 


D*' 


-7r+) = 


= 2/3, 


(5) 


B{D*{2mf - 


D*- 


-7r+) = 


= 1/5. 


(6) 



These branching fractions are assumed throughout. Applying conservation of isospin to 
the spectator decay of Figure 0, it is assumed also that the £)(*)+ £)**o g^j^j j^o _^ 

production rates are equal: 

r(5+ ^ £)(*)+L)**0) = r(s° -> Di*^+D**-). (7) 

This equality is assumed throughout. 

In the Monte Carlo simulation, the three relevant D**^ mass states produce nearly iden- 
tical slow pion 7r~ momentum distributions, resulting in signatures that are virtually indis- 
tinguishable by means of this partial reconstruction technique. For this reason the relative 
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production ratios of Di{j = 1/2)°, Di(2420)°, and D*(2420)° cannot be measured by this 
analysis, but are rather taken from previous experimental results |15;17]- Similarly, it was 
not possible to separate the B D^D** from the B D*^D** modes, and the ratio of 
the branching fractions of these two decays must be assumed. The consequences of both 
assumptions are treated in the discussion of systematic errors. 



C. Partial Reconstruction Kinematics 

In the decays B — > D*^ -D^tt^, the and vr^ are produced nearly 

back-to-back, and the angle a between the reverse direction and vr", shown in Figure |^, 
will be small. For the DgD* and D*D* signals, a ranges between 0° and 30°, with most 
probable values at 11° and 12°, respectively. For D'f^D** signal, a ranges from 0° to 50°, 
with the most probable value at 21°. In contrast to the signal, the background consists of 
uncorrelated and tt^ pairs, for which a will be distributed at random. 
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FIG. 3. Definition of a: the angle between the reverse direction of the measured Df and 
measured vrj. 



It is possible to further constrain a from additional event information that determines 
the allowed D*~ directions. There exist a total of eight unknowns in the decay: the three- 
momentum, the parent D*~ three-momentum, and the two angles governing the B^ direction. 
The B energy is equal to the CLEO beam energy. Requiring conservation of energy and 
momentum in the B^ —>■ DfD*~ and D*~ -D°vr~ decays yields eight constraints, where 
the masses of the B^, Df, D*^ , Z)°, and vr^ are assumed. Solving for the unknowns yields a 
pair of solutions, due to a quadratic ambiguity in the underlying algebra. The procedure 
of this solution follows. 

The D*~ and energies are determined from the measured and energies: 

Ed*- = EBeam — Ejj+, (8) 



Eoo = Ed*--E^^. (9) 

The magnitude of D*~ and momenta follow from their energies Pd*~ = 

^Ejj^_ — M|,*_ and Pdo = ^ Ej^^ — Mj^^ . For the previously- assumed decay, kinematics 
constrain the D*~ to a cone of allowed directions relative to the measured . The radius 
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of this cone, 6i, is shown in Figure ^, and represents the angle between the reverse di- 
rection and inferred D*~ . Using the momentum magnitude, beam energy, and particle 
masses, 6i can be expressed in the lab frame as: 



3280300-005 



Cone of Kinematically 
Allowed D* Momenta 




FIG. 4. Definition of 9i: the angle between the reverse direction and the cone of allowed 
D*^ directions. 



Kinematics also constrain the D*~ to a cone of allowed values about the vr" direction. 
The radius of this cone is 62, the angle between the n~ and inferred D*~, defined in the lab 
frame and shown in Figure ^ 

ML - M2 _ m2 1 
COS 02= ^° ^ , "- +75-5-- 11 
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FIG. 5. Definition of 62'- the angle between the vr^ and cone of allowed D* directions. 



Valid solutions for the D*~ momentum exist at the intersection of the cones defined by 
9i and 6'2, as shown in Figure For the two cones to intersect, the angle between the 
measured and measured vr^ — the angle previously defined as a — must be confined to a 
range bounded by the sum and difference of 6*1 and 6*2: 
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FIG. 6. Valid D* solutions exist where the combined 9i and 62 cones intersect. There are 
generally two solutions, resulting from a quadratic ambiguity in the underlying algebra. 



For a greater than the upper hmit 61 + 62, the smaller cone is entirely outside the larger 
one, preventing their intersection and the existence of a kinematically valid D*~^ solution. 
For a less than the lower limit |6'i — ^2], the smaller cone is completely inside the larger, also 
preventing their intersection. As shown in Figure |^, the lower limit occurs as the smaller 
cone grazes the inside edge of the larger one, where this limit is defined as ag: 



Since only one a and one exist for a particular vr") pair, they are unaffected by the 
quadratic ambiguity in the D*~ solutions. 

In the case of the signal —>■ D^D*~ mode, a is small, ao is as small or smaller, and 
the difference between a and ao is very small. Since the background is relatively isotropic 
in cos a, it is more convenient to work with the cosines of the angles, where it is found 
that the difference cosao — cos a peaks sharply for signal at small values. Signal Monte 
Carlo distributions are shown in Figure § over the range (—0.04, 2.00) for D^D*~ , D*~^D*~, 
and Z)(*)+D**. Two backgrounds are also shown in the figure: the BI3 background, from 
simulated non-signal B meson events, and continuum background, from e~^e~ cc, ss, uu, 
or dd. The three signals display sharp peaking in cos — cos a, where the D^D* and D^;^D** 
peaks are measurably broader than the DgD* . In the D*^D*^ case, this broadness results 
from the random nudge given the by the 7/71° in the D*^ D'^'j/ D^tt^ transition, 
causing the two particles to be not quite so back-to-back. For the D'f^^D**, the broad 
peak comes from the thrust given the D*~ from the unreconstructed tt'^ produced during 
the intermediate D**^ D*^7r+ decay. It should be noted that no sharp peaking occurs 
in either background where the D'^ and vr" are nearly uncorrelated, though some hint of a 
peak is exhibited due to kinematic correlations. 



61-62] <a< 61 + 62. 



(12) 



ao = \6i - 62\. 



(13) 
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FIG. 7. ao is the minimum value a can take for an event where the cones continue to intersect, 
and corresponds to the smaller cone grazing the inside edge of the larger. 



It is seen from Figure ^ that cos ao "~ cos a occasionally drifts below zero. This is due to 
detector resolution effects that distort the quantities used to calculate cos^i and cos 6*2. 



D. Fit of the Data 

A sharp signal peak in the cosoq ~ cos a distribution of the CLEO on-resonance data 
previously described in Section II, superimposed on a relatively flat background, is seen in 
Figure The figure also shows a binned maximum-likelihood fit to this data consisting 
of three components: signal, BB background, and continuum background. The 

signal and BB background components are allowed to float, while the continuum 
level is fixed by scaling the off-resonance background by the on/off- resonance ratio. The 
D{*)D*{*) component is a weighted combination of DgD*, D*D*, and D^*'>D** signals, as the 
signal distribution shapes in cos ao — cos a are too similar for meaningful separation. The 
signal is concentrated in the relatively small region —0.04 < cos ao — cos a < 0.12, where the 
data contains 528 events. Table |T| lists fit results of the three components for this signal 
region. The errors listed are statistical. The subsequent analysis of the relative 
rates and polarizations is confined to the signal region: 

-0.04 < cos ao -cos a < 0.12 (14) 



V. SEPARATION OF DsD% D%D*, AND D^*'^ D** SIGNALS. MEASUREMENT OF 

D*gD* POLARIZATION 

A. Definition of the Two-Dimensional p^^ vs cos6',r Parameter Space 

Once the background levels have been determined, the signal modes may be separated 
from one another. This separation is effected by constructing a two-dimensional parameter 
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FIG. 8. Signal Monte Carlo and background distributions of the partial-reconstruction variable 
cosao - cosQ. Shown are: (a) DsD* Monte Carlo, (b) D*D* Monte Carlo, (c) d''*^ D** Monte 
Carlo, (d) BB background Monte Carlo, (e) continuum data. The signals display a characteristic 
narrow peak, while the backgrounds are relatively broad. 
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FIG. 9. The fitted cosao — cos a data distribution. The fit is broken down into three 
components: D^*^ D*^*'^ signal, BB background, and continuum background. The signal region 
is —0.04 < cosao ~ cos a < 0.12. The continuum background is constrained by scaling the 
off-resonance background by the on/off-resonance ratio. 
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TABLE III. Fit results for the data signal region —0.04 < cosoq ~ cos a < 0.12, containing 
528 events. Errors are statistical. 



Mode 


Number of events 


Signal 


314.0 ±24.0 


BB Background 


138.9 ±5.2 


Continuum 


74.8 (constrained) 



space, where each signal carries a distinctive shape. Two variables are required, of which 
the first is the magnitude of D'^ momentum p^^. The kinematics of the B [)(*)[)*(*) 
decays constrain the relevant Dg momentum range to 1250 MeV/c < pr,^ < 1925 MeV/c. 
The second variable of interest is the cosine of the vr" decay angle as expressed in the D*~ 
frame — cos6'7r — where 6.,^ is a helicity angle, shown in Figure |10[ It is possible to calculate 
cos^^TT from available event information, without reconstructing the D*~: 

2py, 211.00. Mo.p'j,.. ^ ' 




Here the kinematics of the D^D*~ mode have been assumed. The quantities 7^)*, 

Ptt^, and Pd^ are expressed in the lab frame, while E*^^ E^o, and p^o are in the D*~ 
frame, 

p. _ Ml,+Mlo-Ml 

_ Ml, - M^.. + Ml 
2M^. ■ 

In the 5° D^D*~ mode, the D*~ is produced in a (J, J^) = (1, 0) state, and conserva- 
tion of helicity distributes the vr" as cos^ 6'^. Imperfect detector resolution smears the shape. 
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For the case of longitudinally polarized D*~^D*~ from 5° decays, the D*~ is also produced in 
a (1, 0) state. However, the resulting cos^ shape is not centered at the origin, but is rather 
shifted downwards. This shift comes from the missing 7/7?" (where D*~^ Df'y/D^ir^), 
which was not taken into account in the calculation of cos^t^. Nevertheless the original 
cos^ 6',r shape is well preserved, centered at —0.2, and falls over the range (—1.2,0.8). In 
the case of transversely polarized D*^D*^ , the D*'^ is produced in a (1, 1) or (1, —1) state, 
and the resulting vr^ produces a helicity distribution of 1 — cos^Ot,, also centered at —0.2. 
Finally, the three D**^ states each produce the vr" in their unique helicity distributions: the 
Di(2420)° decays as 1 + 3 cos^ 6^, the D*(2460)° decays as 1 - cos^ 6^, and the Di{j = 1/2)° 
decays isotropically. However, blending the three D**^ states according to their production 
ratios in D^D** and D*^D** effectively washes out any characteristic helicity shape. The 
resulting blended distribution is centered at —0.4 and ranges over (—1.4,0.5), because of 
the missing intermediate vr (from D** D*~it), which is not accounted for in calculating 
cos^TT- The limits of cos^^r relevant to this analysis are therefore —1.40 < cos^^ < 1.05. 

The versus cos^^ two-dimensional distributions are shown in Figure ^ for each of 
the four signals {DgD*, Longitudinally Polarized D*D*, Transversely Polarized D*D*, and 
D^*^D**) and two backgrounds {BB and continuum). Because the longitudinally polarized 
and transversely polarized D*D* produce markedly different shapes in this two-dimensional 
distribution, they can be separated into two components. The two-dimensional on-resonance 
CLEO data distribution is shown in Figure |l2l 



B. The Fitted Data 

A two-dimensional binned maximum-likelihood fit is applied to the data. The BB and 
continuum backgrounds, whose levels were determined in the previous one-dimensional fit to 
the cos ao "COS a distribution, are fixed here, and their shapes are parameterized as products 
of Chebyshev polynomials. The two-dimensional Monte Carlo distributions are used for the 
[){*)[)*{*) signals. Four signal components are allowed to float: the number of DsD*, number 
of D*D*, number of D^*^D**, and the relative longitudinal D*D* polarization. Converting 
the likelihood to a x^-like quantity (— 21n£), the resulting fit has a likelihood of 125.4 for 
130 bins with 4 floating parameters. Projections of data and fit along both the p^^ and 
cos^^ axes, broken down into signal and background components, are shown in Figure IT^. 



In Table the number of events resulting from the two-dimensional maximum-likelihood 
fit are reported for each of the D^*^^D*^*^ modes, along with their statistical and systematic 
uncertainties, where the systematics will be discussed in Section |V|. 

The relative longitudinal polarization of the D*^D*~ production is measured to be: 

Fi/F = (50.6 ± 13.9 ±3.6)% (18) 

where the first error is statistical and the second systematic. 

C. The Fitted Data With the D^s^^ D** Component Removed 

Because production oi B [)(*)+[)** j^as not been previously observed, one might 
question its inclusion in the preceding fit. A worthwhile consistency check is to remove the 
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FIG. 11. vs cosSjt distributions for the four signal Monte Carlo and two background 

samples. Shown are: (a) DgD* Monte Carlo, (b) Longitudinally Polarized DID* Monte Carlo, (c) 
Transversely Polarized D*D* Monte Carlo, (d) Di*^ D** Monte Carlo, (e) BB background Monte 
Carlo, (f) continuum data. The box size is proportional to the number of candidates in the bin. 
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FIG. 12. vs cos^tt distributions for the CLEO data. The plot contains 528 events. The 

box size is proportional to the number of candidates in the bin. 

£)(*)+£)** from the set of functions and repeat the two-dimensional fitting procedure. The 
results, without the D'f^^D**, are projected along the pd^ axis and cos6'7r axis in Figure 
The fit projection shape matches the data well, but is shifted systematically upwards by 
about 50 MeV/c. The cos6'7r fit projection shape is decidedly different from the data, since 
the projection is systematically low over the region —1.2 < cos^tt ^ —0.2 and systematically 
high over the region —0.2 < cos^^^ ^ 0.8. The cos^^^r fit projection also displays a pair of 
symmetric peaks that are not refiected in the data. This two-dimensional fit has a likelihood 
of 139.9 for 130 bins with three fioating parameters, and since for this fitting procedure the 
likelihood follows closely the chi"^ behavior, this corresponds to a reduced significance of 3.8 
standard deviations. A study of the data sideband regions off the cos ao — cos a signal peak 
(i.e. where cosoq — cos a > 0.12) reveals an amount of that is consistent with the 

amount of D'f^^D** observed in the signal region —0.04 < cosoq ~ cos a < 0.12. Taking 



TABLE IV. Fitted yield for each DP^ D<*^ mode. The first error is statistical and the second 
is systematic. The D** is the sum of charged and neutral 1)1(2420), D2(2460), and Di{j = 1/2) 
resonances. 



Mode 


Fitted Yield 


S° ^ D+D*- 


92.7 ± 15.3 ±9.5 


D*+D*- 


149.2 ± 30.4 ± 20.9 


B Di*^^D** 


81.6 ± 23.3 ± 15.3 
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FIG. 13. The projections of the two-dimensional CLEO data and fit along the axis (top) 
and cos^TT axis (bottom). The fit contains DgD*, D*D*, Dl D**, and background components. 



20 



the factors all together, we conclude that the data strongly indicate a substantial 
component. 

VI. SYSTEMATIC UNCERTAINTIES 

The single largest uncertainty in the analysis is the 25% uncertainty in the (j)TT~^ 
branching fraction [§: 

B{D^ 07r+) = (3.6 ± 0.9)% (19) 

This uncertainty is displayed separately from the other systematic uncertainties, which are 
listed in Table 0. 

TABLE V. Systematic uncertainties in percent for B — > D^*'^'^ D*''*'' decays and T^/F, the 
longitudinal polarization of Dl'^D*~ . 



Source 


D+D*- 


D*+D*- 






Dg Tracking 


3.0 


3.0 


3.0 


3.0 


TTg Tracking 


5.0 


5.0 


5.0 


5.0 


Total Number of 


1.8 


1.8 


1.8 




BB Mesons 










Fit Normalization 


2.9 


2.9 


2.9 




Monte Carlo Statistics 


1.0 


1.0 


1.0 


1.0 


Continuum Subtraction 


3.7 


3.7 


3.7 


1.8 


BB Background 


1.6 


1.6 


1.6 


1.0 


Subtraction 










Continuum Shape 


1.0 


1.0 


1.0 


1.0 


BB Background 


1.0 


1.0 


1.0 


1.0 


Shape 










Z?i(i = 1/2) :L>i (2420) : 


0.8 


2.9 


4.5 


0.6 


L>|(2460) ratio 










D+D** : D*+D** Ratio 


2.4 


9.8 


14.8 


3.2 


Non-Resonant 


2.2 


4.5 


5.9 


1.2 


Di*^^D*-TT Production 










Total for D^*^^ D*^*^ yield 


9.9 


13.8 


18.5 




B{(I)^K+ K-) 


1.6 


1.6 


1.6 




B{D*- D^TTJ) 


2.0 


2.0 


2.0 




Total systematic uncertainty 


10.2 


14.0 


18.7 


7.3 



A 1% systematic uncertainty in track finding and fitting efficiency is estimated for each 
fast charged track, which for the add linearly to a 3% total. The slow pion tt^ track 
finding and fitting uncertainty is estimated at 5%. The uncertainty in the total number of 
BB meson pairs introduces a systematic error of 1.8%. 
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FIG. 14. Removing the Dl ' D** component from the two-dimensional fit to the CLEO data, 
where data and fit are projected along the pj^+ axis (top) and cos0^ axis (bottom). The fit is 
split into DgD*, D*D*, and background components. The p^^ fit projection shape is shifted 
upwards relative to the data, and the 9.,^ fit projection shape is systematically low over the region 
—1.2 < cos^TT < —0.2 and systematically high over the region —0.2 < cos^jr ^ 0.8. The likelihood 
is reduced by 3.8 standard deviations from the previous fit, which includes a Dg 'D** component. 
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The two-dimensional fit to the data estimates the total amount of signal at 

323.2 events. Since the previous one-dimensional fit to the cos ag— cos a distribution, summa- 
rized in Table |T|, determined the level of D(*)Z)*(*) signal at 314.0±24.0, the two-dimensional 
fit result overestimates the amount of signal by 9.2 events. To test for a systematic bias in 
the two-dimensional fitting procedure, fifty simulated datasets were created and filled with 
[)(*)£)*{*) signal Monte Carlo, BB background, and continuum background according to the 
proportions of Tables |T| and |^ Following the procedure of fixing both backgrounds and 
allowing all four signal components to float, two-dimensional fits to these simulated datasets 
gave fifty estimates of total signal. The difference between the estimate from each 

fit and the number of input events forms a distribution centered at 1.1 with an 

rms of 5.3, consistent with zero and indicative of an unbiased fitting procedure. In the 
case of the fit to the real dataset, the additional 9.2 events differ from the expected total 
by an acceptable 1.7 standard deviations. In order that these events might be accounted 
for, a systematic error of 2.9% is introduced into the overall signal yield. The polarization 
measurement is not affected by this systematic error in the fit normalization. 

Forty thousand signal Monte Carlo events were generated for each of the nine signal 
modes: DfD*~, longitudinally polarized D*^D*~, transversely polarized D*~^D*~, DfD**^ 
(for each of Di{2A20f, D*{2AQ0f, and Di{j = 1/2)°) and D*+D**^ (also for all three Z)**° 
states). To estimate statistical limitations, the signal samples were divided in half and the 
half-samples used to refit the two-dimensional vs cos 9^, data distribution. The resulting 
fits differ from the original by less than 1.0%. 

An uncertainty is introduced by statistical fluctuations in the amount of continuum 
background. Varying the number of continuum background events by one standard deviation 
(cr) affects the overall two-dimensional flt yields by a maximum of 3.7% and the polarization 
by a maximum of 1.8%. The uncertainty from statistical fluctuations in the total number 
of BB background events is anti-correlated with the continuum background. This is the 
result of highly similar background shapes in the one-dimensional flt to the cos — cos a 
data distribution. Refltting the two-dimensional data distribution with these fluctuations 
changes the yields by a maximum of 1.6%, and the polarization by a maximum of 1.0%, 
where the small uncertainty results from the anticorrelation. 

The two-dimensional continuum and background shapes are parameterized as products 
of Chebyshev polynomials. Varying the polynomial coefficients by the parameterization 
errors and refltting the two-dimensional data distribution changes the results by less than 
1.0% for either background. 

In the two-dimensional flt to the data distribution, there is a single component con- 
taining 

signal. The D** label denotes the sum of three L = 1 charm states: 
the Di{j = 1/2), Di(2420), and the D2{2460). Each of these three states has a unique 
mass and width, and produces a different pattern of vr" helicities. In building the D** 
signal component, it is assumed that the D** production rate from B mesons is at a 
Di{j = 1/2) : L)i(2420) : £'^(2460) ratio of 2 : 1 : 6.7, in accordance with the known 
Di(2420) and 1^2(2460) production rates in B ^ D**7r [13] and the preliminary evidence for 
the Di{j = 1/2) []T3|. To understand the systematic bias introduced by this choice of ratios, 
the data was reflt using widely varying ratios of (2 : 1 : 3.3, 2:1: 13.5, 1:1: 6.7, 4:1: 6.7, 
2:2: 6.7, and 2 : 0.5 : 6.7). This caused the D^D* yield to vary by 0.8%, the D*D* yield to 
vary by as much as 2.9%, the D^/^D** yield to vary by 4.5%, and the polarization to vary 
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by 0.6%. 

The D+D** : D*+D** ratio in the D'fW** component has been fixed a priori at 1 : 2 
in the two-dimensional fit. The assumption of this ratio follows from the analogous modes 
B D'fW*, where D^D* : D*D* has been previously measured at 1 : 2, a ratio confirmed 
by this analysis. The pseudoscalar /vector Ds : D* ratio in this spectator decay implies that 
the same ratio should hold for the D^*'>D** case as well. All the decays 5° /}(*)+/)*(*) are 
spectator decays described by a single Feynman diagram and differentiated only by the final 
angular momentum states of the cs [D^ or -D*'*') and cq {D* or D**) quark pairs. To be 
particularly conservative, the DgD** : £)*£)** ratio is allowed to vary between 1 : 1 and 1 : 4. 
Pseudoscalar/vector spin considerations strongly suggest that the ratio be confined between 
these two limits. Varying the D^D** : D*^D** ratio between 1 : 1 and 1 : 4 changes the 
fit results significantly, as the D^D*^ varies by a maximum of 2.4%, the D*~^D*~ by 9.8%, 
the D^*^+D** by 14.8%, and the polarization by 3.2%. These errors are the second largest 
systematic uncertainty, after the — >■ 07r+ branching fraction uncertainty. 

There exists the possibility that significant non-resonant B D'f^^D*n produc- 
tion could contribute to the data sample. The three-body B D^*^^D*7r decay 
peaks nearly as strongly as resonant signal in cosao — cos a. While no measurements 
of the B vr non-resonant production have been made, an analogy can be 

drawn to non-resonant production of -B ^ D*~ttIv. ALEPH has measured the inclu- 
sive branching fraction B D*^7r£^z/ at (1.25 ± 0.25)%, and the product of exclusive 
branching fractions B{B Di(2420)rz/)i3(L>i(2420)° ^ D*7t) = (0.52 ± 0.17)% Jl6[. 
ALEPH has also placed an upper limit on the 1^2(2460) branching fraction at B{B 
D*{2AQ0)£'iy)B{D*{2A60) D*Tr) < 0.39% [0. Although there is no measurement of the 
mode B -Di(j = l/2)£^z/, recent observations at CLEO of the related mode B^ 0**^1^^ 
report B{B L)i(2420)°7r+) ^ 2/3i3(5+ ^ Di(j = l/2)07r+) ^ 2/3i3(5+ ^ D;(2460)07r+) 
Assuming that -B(Di(2420)° ^ D*-ti+) = B{Di{j = 1/2)° ^ D*-n+) = 2/3 and 
i3(L'2(2460)° — * D*"7r+) = 1/5, and assuming that these relative D**n ratios hold in the 
semileptonic case, nearly all of the inclusive B — > D*~t:1~v will be accounted for by resonant 
B — >■ D**l~v. This would leave only a small nonresonant component. Thus a conservative 
upper limit is that non-resonant B D^*^^D*^7r could be as large as 40% of the resonant 
B [)(*)+[)** branching fraction. Three B [)(*)[)** -\- (Non- Resonant) samples were 
created: one that contained 60% pure B with 30% non-resonant B D^D* vr 

and 10% B D*+D*-Tr, one that contained 60% pure B with 10% non- 

resonant B — s> D^D*^7c and 30% non-resonant B D*'^D*^tt, and one that contained 
60% pure B D^;^+D** with 20% non-resonant B D^D* vr and 20% non-resonant 
B D*+D*-n. Refitting the data distribution with these B D^;^+D** + non-resonant 
B D(*)+D*"7r samples changes the results by 2.2% for the B — > DfD*' case, by 4.5% 
for the B D*+D*-, by 5.9% for the B 0^*^+0**, and by 1.2% for the polarization. 
These are the systematic errors listed in Table |V|. Should it be the case that by 60% of 
resonant B [)(*)+[)** branching fraction be non-resonant, the systematic errors increase 
would increase to 3.0% for the B D+D*-, to 6.3% for the B D*+D*-, to 8.1% 
for the B D'f^'^D**, and to 1.8% for the polarization. It should be noted that other 
non-resonant modes, such as B ^ D'^D*~7rn, produce the D'^ in a momentum range that 
is almost entirely below the lower limit of 1250 MeV/c, excluding these modes from this 
analysis. 
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The 1998 PDG values for the and D* branching fractions are — > K^K~ 



(49.1 ±0.8)% and B{D*- D'^n;) = (68.3 ± 1.4)% 111. These introduce systematic errors 
of 1.6% and 2.0%, respectively, into the extraction of the Z}(*)+_D*(*) branching fractions. 

It is assumed in measuring the longitudinal and transverse D*~^D*~ polarizations that 
these final states are independent of one another. In actuality there exists, in the differential 
decay rate, an interference term between the longitudinal and transverse states that is pro- 
portional to the azimuthal angle between the planes of the D*~^ — > -D+7 and D*~ D^7i~ 
decays. This interference vanishes in the integral over the azimuth, and introduces no sys- 
tematic error into the analysis. 

VII. FACTORIZATION AND PREDICTION OF POLARIZATIONS 

The factorization assumption, when expressed in the framework of Heavy Quark Effective 
Theory (HQET) and extrapolating from the form factors measured by the semileptonic B 
decays B D*ii>, allows accurate estimate of the hadronic decay rates for the modes 
B D'^*h, D^*^p, D^;^D^*\ and Additionally, factorization, HQET, 

and the semileptonic decays, predict the relative polarization of the vector-vector hadron 
products for B D*-X decays, such as B ^ D*-p and S° D*+D*- 

We observe a longitudinal polarization in B^ — * D*^D*^ of Tl/T = (50.6 ± 13.9 ± 
3.6)% for = M^,,, where the first error is statistical and the second systematic. The 
observation is consistent with the prediction of (53.5 ± 3.3)% from factorization, HQET, 
and the semileptonic form factor measurements . The same combination also predicts in 
B D*^p a longitudinal polarization of Tl/T = (89.5 ±1.9)% at = M^, which compares 
favorably with the most recent measurement of (87.8 ± 5.3)% ||23|. Finally, predictions are 
also made that at low the longitudinal polarization will be nearly 100%, and at = q^ax 
decreases to 33% 0. Longitudinal polarization as a function of q^ is plotted in Figure |T5| for 



the factorization prediction and compared with the D*~ p and D*~^D*~ measurements. The 
agreement is excellent, confirming the validity of the factorization assumption and HQET 
in extrapolating the semileptonic form factor results for regions of high g^. The polarization 
of the semileptonic decays remains unobserved [p2| . 

Another vector-vector hadronic B decay mode which may further test the factorization 
assumption at high g^ is 5° — *• D*~^D*^ . This decay is Cabibbo-suppressed, and a polariza- 
tion measurement will require higher statistics than those provided by present experiments 
2^. Future experiments will also reduce the errors of the D*~ p and D*^D*~ measurements. 



VIII. SUMMARY AND CONCLUSIONS 

Removing the signal component from the two-dimensional fit reduces the like- 

lihood by 3.8cr, and the resulting projections along both the and cos^^r axes are system- 
atically different from the data as discussed in Section [V^ . Furthermore, a level of 



is observed in the data sideband regions of cos ao — cos a consistent with that seen in the 
signal region. We conclude that the data support first evidence for B — >• D^*^^D** decays. 

From the event yield of Table |rV|, we can calculate the exclusive branching fractions 
5° ^ DfD*-, B° D*+D*-, and 5+ ^ /)(*)+£)**o^ ^j^^^^ ^^le D**^ is the sum of the 
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FIG. 15. Relative fraction of longitudinal polarization in vector-vector B — > D*~X decays as 
a function of q^, where = M'^, and X is a vector meson. Shown are the 1998 measurement of 
B — > D*~ p, and the B^ — > D*'^D*~ polarization measured here for the first time. The shaded 
region represents the prediction using factorization and Heavy Quark Effective Theory, and ex- 
trapolating from the semileptonic B D*ll.v form factor results. The contour is one standard 
deviation [a). 
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^1(2420)°, ^^(2460)°, and Di{j = 1/2)° states 

B{B^ D*~) = (1.10 ± 0.18 ± 0.11 ± 0.28)%, (20) 

B{B° D*+ D*-) = (1.82 ± 0.37 ± 0.25 ± 0.46)%, (21) 

}3{B+ DW+ D**°) = (2.73 ± 0.78 ± 0.51 ± 0.68)%. (22) 

The first error is statistical, the second systematic, and the third the contribution from 
the uncertainty of the Df (pix^ branching fraction. These 5° D^*'>^D*~ branching 
fractions supersede the previous CLEO measurements |l[]. The extraction of the combined 
£){*)+ f)**o branching fraction is contingent on the assumption of Equation (7), where the 
charged-i? decay rate, — > £)(*)+/)**o^ is presumed equal to the neutral--B decay rate, 
50 _^ The extraction also requires some presumption of the individual D** — >■ 

D*^TT rates, shown in Equations (l)-(6). The assumptions follow from conservation of 
isospin in the spectator B decay of Figure It is further assumed that the production rates 
of B~^ and B^ in T(4S') decays are equal for all branching fraction measurements. 

The relative longitudinal D* polarization in B^ — > D*~^D*~ is measured for the first time 

as: 

y^(5° ^ D*+D*~) = (50.6 ± 13.9 ± 3.6)% (23) 

where the first error is statistical and the second systematic. The measurement is consistent 
with the recent factorization prediction of (53.5 ± 3.3)%, confirming the validity of the 
factorization assumption in the domain of relatively high ||21[| . 

ACKNOWLEDGMENTS 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. I.P.J. Shipsey thanks the NYI program of the NSF, 
M. Selen thanks the PFF program of the NSF, A.H. Mahmood thanks the Texas Advanced 
Research Program, M. Selen and H. Yamamoto thank the OJI program of DOE, M. Selen and 
V. Sharma thank the A. P. Sloan Foundation, M. Selen and V. Sharma thank the Research 
Corporation, F. Blanc thanks the Swiss National Science Foundation, and H. Schwarthoff 
and E. von Toerne thank the Alexander von Humboldt Stiftung for support. This work was 
supported by the National Science Foundation, the U.S. Department of Energy, and the 
Natural Sciences and Engineering Research Council of Canada. 



27 



REFERENCES 



CLEO Collaboration, D. Gibaut et al, Pliys. Rev. D 53, 4734 (1996). 

J.L. Rosner, Phys. Rev. D 41, 3732 (1990). 

A. Falk and M. Luke, Phys. Lett. B 292, 119 (1992). 

M. Neubert, Phys. Lett. B 264, 455 (1991). 

V. Rieckert, Phys. Rev. D 47, 3053 (1993). 

G. Kramer, T. Mannel, and W.F. Palmer, Z. Phys. C 55, 497 (1992). 

ARGUS Collaboration, H. Albrecht et al, Z. Phys. C 48, 543 (1990). 

CLEO Collaboration, M. Artuso et al, Phys. Lett. B 378, 364 (1996). 

CLEO Collaboration, G. Brandenberg et al, Phys. Rev. Lett. 80, 2762 (1998). 

CLEO Collaboration, Y. Kubota et al, Nucl. Instrum. Methods Phys. Res. A 320, 66 

(1992). 

R. Brun et al, CERN Report No. CERN-DD/EE/84-1, 1987 (unpublished). 
Particle Data Group, C. Caso et al, Eur. Phys. J. C 3, 1 (1998). 
CLEO Collaboration, P. Avery et al, Phys. Lett. B 331, 236 (1994). 
E687 Collaboration, S. Frabetti et al, Phys. Rev. Lett. 72, 324 (1994). 
CLEO Collaboration, S. Anderson et al. Report No. CLEO CONE 99-6, XIX Inter- 
national Symposium on Lepton and Photon Interactions at High Energies, Stanford, 
California, 1999 ( |hep-th/9908009|) . 



ALEPH Collaboration, D. Buskuhc et al, Z. Phys. C 73, 601 (1997). 

CLEO Collaboration, T. Bergfeld et al, Phys. Lett. B 340, 194 (1994). 

T. Browder, K. Honscheid, and D. Pedrini, in Annual Review of Nuclear and Particle 

Science, 46, 395 (1996). 

M. Neubert, Phys. Rep. 245, 259 (1994). 

Z. Ligeti, Y. Nir, and M. Neubert, Phys. Rev. D 49, 1302 (1994). 

J. D. Richman, in Probing the Standard Model of Particle Interactions, edited by 

R. Gupta, A. Morel, E. de Rafael, and F. David (Elsevier, Amsterdam, 1999), p. 640. 

CLEO Collaboration, J. Duboscq et al, Phys. Rev. Lett. 76, 3898 (1996). 

CLEO Collaboration, G. Bonvicini et al, CLEO CONE 98-23, ICHEP98 852 (1998). 

CLEO Collaboration, M. Artuso et al, Phys. Rev. Lett. 82, 3020 (1999). 



28 



